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OPTIMAL ENERGY DECAY IN A ONE-DIMENSIONAL 
COUPLED WAVE-HEAT SYSTEM 


CHARLES BATTY, LASSI PAUNONEN, AND DAVID SEIFERT 

Abstract. We study a simple one-dimensional coupled wave-heat sys¬ 
tem and obtain a sharp estimate for the rate of energy decay of classi¬ 
cal solutions. Our approach is based on the asymptotic theory of Co¬ 
semigroups and in particular on a result due to Borichev and Tomilov 
m, which reduces the problem of estimating the rate of energy decay 
to finding a growth bound for the resolvent of the semigroup generator. 
This technique not only leads to an optimal result, it is also simpler 
than the methods used by other authors in similar situations. 


1. Introduction 

We study the following one-dimensional coupled wave-heat system: 



^G (-1,0), t>0 


^G (0,1), t>0 

u^{—l,t) = 0, w{l,t)=0, 

t > 0 

ut{0, t) = w{0, t), u^{0, t) = w^{0, t), 

t > 0 

m(^,0) = m(O, ut{C,0) =v{^), 

^G(-1,0), 

w{C,0) = w{C), 

^G (0,1), 


where the initial data satisfy tt G 1,0), v G L^(—1,0) and tc G L^(0,1). 

Models of this type have received considerable attention in recent years, es¬ 
pecially in the higher-dimensional setting. The main motivation comes from 
the study of so-called fluid-structure models, in which the non-linear elas¬ 
ticity equation is coupled with the Navier-Stokes equations; see for instance 
the survey article [3]. The above model can be viewed as a simple linearised 
version of such a model which however preserves the feature of a hyperbolic 
equation being coupled with a parabolic equation. 

The primary aim in such models is to obtain quantitative estimates for 
the rate of energy decay of a given solution. Given a vector x = (u, v, w) 
of initial data as above, the energy of the solution corresponding to these 
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initial data is defined as 
1 

Exit) = 2 J ^ + \M^:t)\^ + \wi^,t)fd^, t > 0, 

where all functions have been extended by zero in ^ to the interval (—1,1). 
A simple calculation shows that 

Exit) = -[ \wd^,t)\‘^ d(, t>0, 

Jo 

provided the solution is sufficiently regular, and hence the energy of any 
such solution is non-increasing in time. The aim of this paper is to establish 
an optimal estimate on the rate at which this decay occurs. 

A higher-dimensional version of the above problem was studied in m, 
showing that Exit) = 0(t as t ^ oo, and this rate was improved in 
[12] to for any <5 > 0. A rate of is established for a closely related 

model in [7] , where the authors also comment on the optimality of the decay 
rates for the higher-dimensional systems; see [T] Rem. 1.3]. Further related 
results may be found in 13 El da US). As the main result of this paper we 
show that the energy of any classical solution to dni) satisfies 

Exit) = oit~^^, t ^ oo. 

Moreover, we prove that this rate is optimal. A similar result was obtained 
in m for a similar one-dimensional problem in which the wave part satishes 
a Dirichlet boundary condition rather than a Neumann boundary condition 
at ^ = —1. Our method is based on the semigroup approach used in all of 
the aforementioned works; see also El Si- However, rather than estimating 
the norm of the semigroup directly we follow [7] in using a recent result 
in the abstract theory of Co-semigroups obtained by Borichev and Tomilov 
in [To], which makes it possible to deduce sharp rates of energy decay from 
appropriate growth bounds on the norm of the resolvent of the generator; see 
also urn- This method not only leads to an optimal decay estimate for the 
problem under consideration, it also leads to a relatively concise argument. 
In particular, the present method appears to be significantly simpler than 
the method used in the Dirichlet case in [16], which is based on the theory 
of Riesz spectral operators and a very detailed spectral analysis. 

The paper is set out as follows. In Section [2] we show that our problem is 
well-posed in the sense of Co-semigroups and we describe the spectrum of 
the inhnitesimal generator. In Section E| we turn to resolvent bounds. We 
first establish in Theorem 13.11 an upper bound for the norm of the resolvent 
operator along the imaginary axis. This is the most technical part of the 
paper. We then go on to show in Theorem 13.41 that the upper bound is 
optimal. In Section 0] we apply the Borichev-Tomilov result to deduce from 
these resolvent bounds an optimal estimate for the rate of energy decay of 
smooth solutions. Finally, in Section [5] we conclude with some comments 
on the Dirichlet case addressed in |16j , explaining how the present approach 
can be adapted to that setting without difficulty. 

The notation used is standard throughout. In particular, given a closed 
operator A on a Banach space X, which will always be assumed to be 
complex, we denote its domain by D(A), its kernel by Ker(A) and its range 
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by Ran(A). The spectrum of A is denoted by and its resolvent set by 

p(A). Given A E p(A), we write R(A, A) for the resolvent operator (A — 

Given two functions f,g : (0,oo) —>■ M+, we write f{t) = 0{g{t)), t —>■ oo, 
to indicate that f{t) < Cg{t) for some constant C > 0 and all sufficiently 
large t > 0. If g{t) > 0 for all t > 0, we write f{t) = o{g{t)), t —>■ oo, 
if f{t)/g{t) —0 as t —)• oo. The case of real-valued functions defined on 
the whole of M is treated in an analogous way. For real-valued quantities p 
and q, we occasionally use the notation p ^ q to indicate that p < Cq for 
some constant G > 0 which is independent of all the parameters that are 
free to vary in a given situation. We let C_ denote the open left half-plane 
{A E C : ReA < 0}. 

2. WeLL-POSEDNESS and PROPERTIES OF THE SEMIGROUP 

We begin by recasting dni) as an abstract Gauchy problem. Gonsider 
the two Hilbert spaces X = x L^(—1,0) x L^(0,1) and Y = 

X 1), both with their natural norms. The op¬ 

erator A defined by Ax = {v, u", w") for x = (u, u, w) in the domain 

D{A) = {{u,v,w) E Y : u'{—l) = u;(l) = 0, u(0) = tc(0), u'(0) = tc'(O)} 

is closed and densely defined. If we let z{t) = (u(-, t), ut(-, t), !/;(•, t)) for 
t > 0, then the coupled wave-heat system a can be rewritten in the form 


z'{t) = Az{t), t > 0, 
2;(0) = X, 


where x G X. We show below, in Theorem l2.ll that A generates a uniformly 
bounded Co-semigroup {T{t))t>o on X. In particular, the unique solution 
of (12.ip is given by z{t) = T{t)x for t > 0. This solution in general satisfies 
only in the so-called mild sense, and it is a solution in the classical 
sense precisely when x E D{A)\ see for instance [T] for details on the theory 
of Co-semigroups. It will be for such classical solutions that we eventually, 
in Section 01 obtain a uniform rate of energy decay in the sense described in 
Section [H 


Theorem 2.1. The operator A generates a uniformly bounded Co-semigroup 
{T{t))t>o on X. Moreover, the spectrum cr{A) of A consists of isolated 
eigenvalues and satisfies cr{A) n iM = {0} and, in fact, 

cr{A) = {a E C : \/Acosh(A) cosh(\/A) -|-sinh(A) sinh(\/A) = O}. 

We prove this theorem in a sequence of smaller results. 

Proposition 2.2. The spectrum (t{A) of A consists of isolated eigenvalues 
and satisfies (t{A) n = {0}. Moreover, Ker(H) = span{(l, 0,0)}. 

Proof. By the Rellich-Kondrachov theorem Y embeds compactly into X. 
Since the resolvent operator R{X,A), for any A E p{A), maps X isomor- 
phically onto D{A) endowed with the graph norm and this space embeds 
continuously into Y, A has compact resolvent. In particular, the spectrum 
of R{\, A) consists only of eigenvalues of finite multiplicity whose only pos¬ 
sible accumulation point is the origin. By the spectral mapping theorem for 
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the point spectrum, cr^A) consists only of eigenvalues of finite multiplicity 
whose only possible accumulation point is at infinity. 

For elements x = {u, v, w) and y = (/, g, h) of X, let 

(x I y) := {u', f)L 2 + {v, g)L 2 + {w, h)L 2 . 

Suppose that s G M with s 7 ^ 0 and that x = {u,v,w) G Ker(is — A), i.e. 
X G D{A) and {is — A)x = 0. A routine calculation using integration by 
parts gives 

0 = Re {{is — A)x | x) = Re {—Ax \ x) = ||rc^||^ 2 . 

Since w G H^{Q,1) and rc(l) = 0 this implies that re = 0. The equation 
{is — A)x = 0 now reduces to 

( 2 . 2 a) u"{0 = -M0, eG(-l, 0 ), 

( 2 . 2 b) v{i) = isu{£,), |g(- 1 , 0 ), 

( 2 . 2 c) «'(-!) = u'{ 0 ) = v{ 0 ) = 0 . 

The general solution of (I2.2ap with the boundary condition u^(—1) = 0 takes 
the form 

u (0 = o(s)cos(s(C + 1 )), Cg[- 1 , 0 ], 

where a{s) G C, and the boundary conditions at ^ = 0 give 
a{s) sin(s) = is a{s) cos(s) = 0 . 

Since s ^ 0 these two equations imply that a(s) = 0. Thus the solution of 
(12. 2 p is u = x = 0, and hence x = 0. It follows that Ker(zs — A) = {0} for 
s 7 ^ 0, so that o'p{A) n C {0}. 

It remains to consider the case s = 0. Certainly (c, 0,0) G Ker(A) for any 
c G C. On the other hand, if x = {u,v,w) G Ker(A), then u; = 0 as above 
and 

x"( 0 = 0 , ^€(- 1 , 0 ), 

x( 0 = 0 , ^€(- 1 , 0 ), 

u'(—1) = x(0) = u'(0) = 0, 

whose solutions are of the form u = c, for some constant c G C, and x = 0. 
Thus 0 G o'p{A) and Ker(A) = span{(l,0,0)}, as required. □ 

Proposition 2.3. The eigenvalues of A are precisely those points A G C 
which satisfy 

(2.3) \/Acosh(A) cosh(\/A) + sinh(A) sinh(\/A) = 0. 

Proof. Suppose that A 7 ^ 0 and let x = {u, x, w) G Ker(A — A). Then 
(2.4a) x"(0 = AV0, Cg(-1,0), 

(2.4b) ^;(e) = An(e), ^G(-1,0), 

(2.4c) w"{f,) = Xw{f,), (, G (0,1), 

(2.4d) u'{—l) = tx(l) = 0, x(0) = tx(0), x'(0) = rx'(O). 

The general solution of (|2.4al) with the boundary condition u'(—1) = 0 is 
(2.5) n(0 = a(A)cosh(A(^ +1)), ^g[-1,0]. 
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where a(A) G C. Moreover, the solution of ()2.4cl) with the boundary condi¬ 
tion rc(l) = 0 is given by 

(2.6) w{^) = b{X)smh{VX{^-l)), ^g[0, 1], 

where 6 (A) G C. Using (|2.5I) and (|2.6I) . the coupling condition w{0) = u(0) = 
An(0) becomes 

Aa(A) cosh(A) -|- 6 (A) sinh(\/A) = 0. 

Similarly, the condition w'{0) = tt^(O) can be written in the form 
Aa(A) sinh(A) — \/A 6 (A) cos h(VA) = 0. 

Thus Ker(A — A) = {0} if and only if the equation 

/Acosh(A) sinh(VA) \ /a(A)\ _ ^ 

\Asinh(A) — \/Acosh(\/A)/ 

has only the trivial solution a(A) = 6 (A) = 0. In other words, A G (t{A) 
if and only if the matrix on the left-hand side has zero determinant. For 
A 7 ^ 0, this is equivalent to (|2.3I) . □ 

Proposition 2.4. The space X splits as a topological direct sum 

X = Ran( 74 ) 0 Ker(A), 

and on the closed subspace Ran(^) the norm of X is equivalent to the norm 
III • III given by 

|||(r,u,u;)||| = (||u '||^2 + \\v\\\2 + Ik||i 2 )^^^ . 

Proof. We begin by showing that Ran(^) = Ker((/>), where f : X ^ C is 
the bounded linear functional that maps x = {u, v,w) G X to 

(2.7) fix) = r(0) 0 j ^ u(0 (1 “ dC 

This shows in particular that Ran(^) is closed and that the projection x i-A 
((^(x),0,0) onto Ker(j4) along Ran( 74 ) is bounded. 

Suppose that y = {f,g,h) G Ran( 74 ). Then there exists x = {u,v,w) G 
D{A) such that Ax = y, which leads to the equations 

( 2 . 8 a) Cg(- 1 , 0 ), 

( 2 . 8 b) v{0 = f{0, ^€(- 1 , 0 ), 

( 2 . 8 c) w"i^) = h{0, CG( 0 , 1 ), 

(2.8d) u'{-l) = t(;(l) = 0, u(0) = w{0), 'u'(O) = u)'(O). 

A direct computation shows that the solutions of (I2.8ap - (l2.8cl) . subject to 
the boundary conditions u'{—l) = 0 and rc(l) = 0 , are of the form 

u{^)=a + j {C-r)g{r)dr, 


HO 


fiO, 

fe(l -0 + ^ (r-OKOdr, 








6 


CHARLES BATTY, LASSI PAUNONEN, AND DAVID SEIFERT 


where a,b £ C are constants. It is now easy to verify that the coupling 
conditions u(0) = r(;(0) and tt'(O) = can both be satisfied only if 

0(y) = 0, which shows that Ran(^) C Ker((/)). 

Conversely, if y = {f,g,h) E X is such = 0, then we can define 
X = {u, v,w) gY by 

= j iC- r)gir) dr, 

v{C) = /(O, 

wiO = (1 - 0 (/(O) - i: rh(r) dr'^ + J (r — ^)h{r) dr. 

A direct computation shows that x E D{A) and that u,v,w satisfy (|2.8I) . 
Thus Ax = y, and hence Kev{4>) C Ran(A). 

By Proposition 12.21 Ker(^) = span{(l,0,0)}. Together with the above 
characterisation of Ran(^) this shows that X = Ran(A) © Ker(A). It re¬ 
mains to show that the norms ||| • ||| and H-H are equivalent on Ran(A). Let 
X = {u,v,w) E Ran(X). Then 

fO 

\\u\\i 2 < sup |u(0| < |u(0)| + sup / |u'(r)|dr. 

Since x E Ker(</>), 

|n(0)| < / + 

and hence ||u ||^2 < ||u^||l 2 + ||'y||L 2 + ||r(;||^ 2 . It follows that |||x||| < ||x||. 
Since ||x|| < |||x||| trivially for all x £ X, the proof is complete. □ 

Proof of Theorem 12.11 We show first that A generates a Co-semigroup on 
X. Let X = {u,v,w) £ D{A). A straightforward computation using the 
definition of D{A) and the Cauchy-Schwarz inequality shows that 

Re((A — I)x, x) < ||u||i 2 ||u ||/^2 — IIre'll ^2 — ||n||^i — 11^11^2 — ||re||^ 2 , 

and it follows from the scalar inequality ab < + 6^) that A — / is 

dissipative. Since 1 E p{A — I) by Propositions 12.11 and 12.31 the Lumer- 
Phillips theorem shows that A — I generates a Co-semigroup of contractions 
on X, and hence A generates a Co-semigroup {T(t))t>o. H remains to show 
that {T{t))t>o is uniformly bounded. 

By Proposition 12.41 X = Xo © Xi topologically, where Xo = Ran(A) and 
Xi = Ker(A) are closed (r(t))t>o-invariant subspaces of X. Moreover, Xo is 
a Hilbert space with the equivalent norm ||| ■ ||| induced by the inner product 

{x\y) = («', /')i 2 + {v, g)L 2 + {w, h)L 2 . 

where x = {u, v, w) and y = (/, g, h). For j = 0,1, let Aj denote the restric¬ 
tion of A to Xj, so that D{Aj) = D{A) n Xj. Furthermore, let {Tj{t))t>o 
denote the Co-semigroup obtained by restricting {T{t))t>o to Xj, noting that 
Aj is the infinitesimal generator of {Tj{t))t>o. Given x = {u,v,w) £ D{Aq), 
then similarly as in the proof of Proposition 12.21 
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Thus Aq is dissipative. Note also that the splitting of the space coincides 
with the spectral decomposition for the eigenvalue 0. In particular, a{Ao) = 
(t(A)\{ 0}, and hence 1 G p{Aq). It follows from the Lumer-Phillips theorem 
that (To(^))t>o is a contraction semigroup. Moreover, for x ^ Xi and t > 0, 
Ti{t)x = X. Hence, given x £ X and x = xq + xi with xq £ Xq, xi £ Xi, 

ll^(^)a;|| < ll|7b(i)a;o||| + ||a;i|| < |||xo||| + lki|| < ||x||, t > 0, 

where the implicit constants are independent of t. It follows that {T{t))t>o 
is uniformly bounded, as required. □ 

Remark 2.5. The above proof shows that the restriction Aq of A to Ran(H) 
satisfies (t{Aq) = it(j 4)\{0}. In particular, (t(j4o) C C_. 


3. Resolvent estimates 


In this section we study the the behaviour of the resolvent operator 
R{is,A) as |s| —>■ oo, where A is the generator of the semigroup {T(t))t>o 
studied in Section [21 First, in Section IXTl we obtain an upper bound on 
the growth of ||R(is,H)|| as |s| —)• oo, and in Sectionwe show that this 
estimate is optimal. These results will allow us to deduce sharp estimates 
on the rate of energy decay in Section U] below. 

3.1. An upper bound. The main result of this section is the following 
asymptotic upper bound on the operator norm of the resolvent operator. 

Theorem 3.1. We have ||i?(is, A)|| = 0(|s|^/^) as |s| —)• oo. 

We begin with two technical estimates. 

Lemma 3.2. There exists c > 0 such that 


is cos(s) cosh(\/is) + i sin(s) sinh(\/is) 


> cexp 


V2 


|s| > 2. 


Proof. Let s G M with |s| > 2. By splitting y/is into real and imaginary 
parts and writing the hyperbolic functions in terms of exponentials, it is 
easy to see that 

|S|V2' 


2 exp — 


V2 


is cos(s) cosh(-\As) + i sin(s) sinh(\/is) 


> |^/^scos(s) ± zsin(s)| — exp ( — \/2|s|^'^^)|\/iscos(s) =F isin(s)|. 
Note that 

2 |s| 

|^/^scos(s) ± i sin(s)| =-^ cos^(s) + f —^ cos(s) ± sin(s)J . 

Thus if |s|^/^|cos(s)| > 1/2, then 

(3.1) |\/is cos(s) ± zsin(s)| > — 

Suppose, on the other hand, that |s|^/^|cos(s)| < 1/2. Since |s| > 2, it 
follows that |cos(s)| < l/(2\/5), and hence |sin(s)| > l/\/2. Thus 

|S|V2 


V2 


■ cos(s) ± sin(s) 


|s|V2 

> |sm(s)| - -^|cos(s)| > 


2^2’ 















CHARLES BATTY, LASSI PAUNONEN, AND DAVID SEIFERT 


SO dsn holds in this case as well. Since 

exp ( - |\/is cos(s) ^ i sin(s )|<^ 

whenever |s| > 2, the result follows. 


3 ^ 2 V 2 


□ 


Lemma 3.3. There exists a constant C > 0 such that, for all f £ 1,0), 

g E L^(—1,0), s E M \ {0} and f E [—1,0], 

sin(s(^ — r)) (zs/(r) + g{r)^ dr 


'-1 


J cos{s{f, — r))(^isf{r) + g{r)^ dr 


^ C'II/IIj^i + 11511^2, 

<C\\f\\m + \\gh2. 


Proof. Using integration by parts it is easy to show that 


L 


; /■? 

sm{s{f, — r))sf{r) dr = / (l — cos(s(^ — r)))/'(r) dr 

1 J~i 

+ (1 - cos(s(C + l)))/(-l). 


Hi- 


Since the evaluation functional / i-A /(—I) is continuous on 1,0), it 

follows that 

/■? 

/ sin(s(^ — r))s/(r) dr 
An analogous argument shows that 

cos(s(^ — r))s/(r) dr 


'-1 


< 


Hi) 


and the result follows by an application of the Cauchy-Schwarz inequality 
to the terms involving g. □ 


Proof of Theorem 13.11 For s E M with |s| > 2, let y = {f,g,h) E X and 
define x = {u,v,w) E D{A) by x = R{is,A)y. Then 

Ikll = {\\u\\hi + ||ii|li2 + \\w\\l2) < \\u\\l2 + ||n'||i2 + \\v\\l2 + \\w\\l2 

and since {is — A)x = y implies that v = isu — /, this becomes 
(3-2) ||x|| < ||su||i2 + ||n'||£,2 + ||u;||i2 + 11/11^2, |s| > 2. 

Thus the result will follow once we have established the following estimates: 
(3-3a) IIs'wIIlz < |s|^'^^||/||//i + 1-51^^^115111,2 + ||/i||^2, 

(3-3b) II'w1Il2 < |'S|^'^^||/||hi + I'5|^'^^||5||l2 + II^IIl 2) 

(3-3c) II'mi||l 2 ^ II/IIhi + II 5 ||l 2 + II^IIl 2 

for all s E M with |s| > 2. Indeed, by (|3.2|) these estimates imply that 
111^11 ~ I lid II for |s| > 2. Since 5 E A was arbitrary, it then follows that 
||dZ(is, A)|| = 0(|s|^/^) as |s| —>■ 00, as required. 

Let s E M be such that |s| > 2. We begin by deriving formulas for 
the components u and w of the vector x = {u,v,w). Note first that the 
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equation (is — A)x = y is equivalent to the following system of boundary 
value problems: 

(3.4a) = -s'^uiC) - isfiO - £/(0> C G (-1,0), 

(3.4b) v{^) = isu{^) - f{^), ^€(-1,0), 

(3.4c) w"(0 = isw{^) - h{^), ^ E (0,1), 

(3.4d) ^^(—1) = ^(1) = 0) 'f^(O) = '*^(0)) u'{0) = w'(0). 


Let 

Us{0 = -[ sm{s{C-r)){isf{r)+g{r))dr, ^e[-1,0], 
■s J-l 

and note that 


U'siO = j ^ cos(s(C - r)) {isf{r) + g{r)) dr, 


[- 1 , 0 ]. 


The general solution of the differential equation (|3.4al) with the boundary 
condition u'{—l) =0 can be written as 


(3.5) n(0 = a(s)cos(s(C +1)) - [/^(O, Ce[-1,0], 


where a(s) E C. In particular, 

(3.6) n'(0 =-sa(s)sin(s(^ +1)) - f7'(0, ?E[-1,0]. 


Furthermore, let 
(3.7) WsiC) 


Sin. 


IS 


h{Vis{r — ^))h{r)dr, 


[ 0 , 1 ], 


noting that 

(3.8) 1T'(0 = i: cosh(\/is(r — ^))h{r)dr, ^ E [0,1]. 

The general solution of the differential equation (j3.4cD with the boundary 
condition in(l) = 0 is given by 

(3.9) w{^) =-b{s)sinh{Vis{l - ^)) +Ws{^), ^E[0, 1], 

where b{s) E C. In particular, 

w'{^) = Vis b{s) cosh(\/is(l - 0) + ^ [0, !]• 

It remains to determine the constants a{s) and b{s). Using (13.511 and (|3.9p . 
the coupling condition r(;(0) = n(0) = zsu(O) — /(O) is easily seen to be 
equivalent to 

isa{s) cos(s) + b{s) sinh(\/is) = /(O) + isUs{0) + 114(0), 
and similarly the condition w'{0) = u'{0) is equivalent to 

sa{s) sin(s) + Vis b{s) cos h(Vil) = -u'M - iT'(o). 

These two equations can be written in the form 
/iscos(s) sinh(Vis) \ /a{s) 

\ssin(s) -\/is cosh (Vis)/ 



(3.10) 


f{0) + isUsi0) + Ws{0) 

. -U'M - lu'(o) 







10 


CHARLES BATTY, LASSI PAUNONEN, AND DAVID SEIFERT 


Let 


Then 


M{s) = 


'iscos{s) sinh(\/is) \ 
^ssin(s) y/is cosh.{y/is)) 


(3.11) det M(s) = cos(s) cosh(\/is) — s sin(s) sinh(\/is), 

and Theorem 12.11 implies that detM(s) ^ 0 for all s G M with s 7^ 0. In 
particular, M{s) is invertible for |s| > 2. Solving (I3.10p finally gives 


(3.12) 


(3.13) 


+ =L^(C(o)+ »''(»)), 


b{s) = - 


det M {s 
s sin(s) 


detM(s) 


{fi0) + isUsi0) + Ws{0)) 


is cos(s) 


{u'M + w'{o)). 


det M(s) 

We now establish the estimates (|3.3al) and (I3.3bl) . By (|3.5p and (|3.6p . 

(3.14) \su{^)\ < |sa(s)| + |sC/s(OI and |m'( 0I < |'So('S)| + \U's{0\ 


for all ^ G [—1,0], and Lemma [331 gives 

(3.15) sup |s[/s(OI, sup |L'^(0I ^ II/IIhi + II5||l2- 

«G[-i,o] ?e[-i,o] 

It therefore remains to estimate |sa(s)| for jsj > 2. A simple estimate based 
on Lemma 13.21 shows that 


(3.16) 


s cosh(r\/is) 


s sinh(r\/Is) 

det M(s) 

rO cIaJ-U 

det M(s) 


<1 


for jrj < 1 and jsj > 2. Using these together with (|3.7p . (|3.8I) and (j3.15p in 
(|3.12p gives 


|sa(s)| < 


scosh(\/is) / s sinh(\/is) / n 


det M{s 

^ s sinh(\/is(l — r)) 


+ 


det M(s) 


det M(s) 
h{r) dr 


< l.s|V2 


m 


+ l'S|^'^^llfl'llL2 + \\h\\l2 


for jsj > 2. Thus (I3.3al] and and (l3.3bD follow from (|3.14D and (|3.15D . 

We now turn to proving (I3.3cp . By p3.9p , (13.111) and (I3.13p , we have that 

w{C) = ^ (sin(s) (/(O) + isUs{0)) + z cos(s)[/'(0)) 

+ (sinh(Vis(l - O)^s(O) - sinh(v^)IUs(0) 

det^M{s) “ O)^s(O) + \/iscosh(Vis)IUs(0) 
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for all ^ E [0)1]- Let iPi(^), W 2 {^), W 3 {^), respectively, denote the three 
terms on the right-hand side of this equation. We estimate wi, W 2 , W 3 in 
turn. From (j3.15l) and (13.161) it is clear that 

l'iei(OI ^ ll/lliLi + IbllL^j C £ [O,!]- 

Turning to W 2 , it follows from an elementary calculation using (13.71) that 
2 -\As(^sinh(\/is(l — ^))lFs( 0 ) — sinh(\/is)lFs(^)^ 

= / cosh.{V^(l — ^ — r))h{r) dr — / cosh(-\/is(l — ^ + r))/i(r) dr 

Jo Jo 

1 

cosh(\/is(l + ^ — r))h{r) dr 
for all ^ E [0,1]. Hence a simple estimate using (j3.16l) gives 

|s|'/> 2 ( 6 I <l|/i|lL^, Ce[o,i]. 

An analogous argument shows that 

^ 3(01 ^ I|L||i, 2 , Cg[o, 1 ]. 

Combining these estimates gives (I3.3cp . thus completing the proof. □ 

3.2. Optimality of the upper bound. The following result shows that 
the spectrum (t{A) of A contains sequences of eigenvalues whose real parts 
tend to ±00 and whose imaginary parts approach zero at a rate which shows 
that the resolvent estimate in Theorem 13.11 cannot be improved. 

Theorem 3.4. There exist sequences (A^) in cr{A) and a constant C > 0 
such that Im ~ ±n 7 r as n —)• 00 and 

(3.17) < ReA^ < 0, n > 0. 

I ImAn 

In particular, 

limsup |s|“^'^^||i?(is, A)|| > 0 . 

|s |—>-00 

Proof. Consider the meromorphic functions F and G defined by F(X) = 
coth(A) and G{\) = tanh(\/A)/\/A, where the square root is defined with a 
branch cut along the negative real axis. By Theorem 12.11 the roots of the 
function F + G are eigenvalues of A. We use Rouche’s theorem to determine 
the approximate location of such roots; see for instance [ 11 ] and [T 6 | for 
related arguments. 

The function F is vri-periodic, with simple poles at A = mri and simple 
zeros at A = Xn, where Xn = {n+^)TTi, n E Z. Since F'{Xn) = 1 for all n E Z, 
a simple argument using a Taylor expansion shows that |F(A)| > ^|A — A^l 
provided |A — A„| is sufficiently small. For n E Z, let rn = \n + and 

= {A E C : |A — A,i| < 2r„}. Then for all sufficiently large \n\ we have 
that |F(A)| > for A E and both F and G are holomorphic in a region 
containing the closure of Simple estimates show that 

sup |G(A)| ~ n —>■ ± 00 , 
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and hence |G(A)| < |-F(A)| for A E dVtn with |n| sufficiently large. By 
Rouche’s theorem, F and F + G have the same number of zeros inside 
for n E Z as above, namely one. Theorem 12.II implies that (t(A) C C_U{0}, 
so (j3.17l) follows. The final statement is a simple consequence of the fact 
that ||R(is,^)|| > dist(zs, < 7 ( 74 ))“^ for all s 7 ^ 0. □ 

4. Energy decay 

Theorem 32] below gives a quantified estimate for the rate of energy decay 
of classical solutions to (ini). The proof relies on the following abstract 
result; see [ini Theorem 2.4]. 

Theorem 4.1. Let {T{t))t>o be a uniformly bounded Co-semigroup on a 
Hilbert spaee X. Let A be the generator of {T{t))t>o and suppose that 
cr{A) n = 0. Then, for any eonstant a > 0, the following conditions 
are equivalent: 

(i) ||i?(is,^)|| = 0(|s|") as jsj — 00 ; 

(ii) ||T(t)^“^|| = as t ^ 00 ; 

(hi) ||T(t)x|| = as t —)• 00 for all x E D{A). 

The following result is the main result of this section. 

Theorem 4.2. If x ^ D{A), then Ex{t) = o{t~'^) as t ^ 00 . 

Proof. Using the notation introduced in the proof of Theorem 12.11 

= ^iro(t)xor, t>o, 

where x = xq xi with xq E Xq and xi E Xi. For y G Xq, it follows from 
the equivalence of the norms || • || and ||| • ||| on Xq that 

|||R(is,^o)y||| < \\R{is,A)y\\ < ||R(is, A)|||||y|||, s / 0, 
and hence by Theorem 13.11 

||R(is,^o)|| = |s|-;^oo, 

where the operator norm is taken with respect to the ||| • |||-norm on Xq. 
Since < 7 (^ 0 ) H zM = 0 by Remark 12.51 it follows from Theorem 14.11 that 
|||ro(t)xo||| = as t —>■ 00 , and the result follows. □ 

Remark 4.3. (a) The rate t~^ in Theorem l4.2l is optimal in the sense that, 
given any positive function r satisfying r{t) = o{t~‘^) as t —>■ 00 , there 
exists X E D{A) such that Ex{t) / o{r{t)) as t —)• 00 . This follows from 
Theorem 13.41 and the uniform boundedness principle together with [9l 
Proposition 1.3]; see also [U Theorem 4.4.14]. 

(b) If we denote by Aq the restriction of A to Ran(^) as in the proof of 
Theorem EH then Theorem 13.41 implies that the spectral bound of Aq 
satishes s(Ao) = 0. Since (To(t))t>o is uniformly bounded, its growth 
bound is equal to a;o(Ro) = 0. Suppose there exists a positive function 
r such that r(t) —>■ 0 as t —)> 00 and for which Ex{t) = 0{r{t)) as 
t —)• 00 for all X E X. Then by the uniform boundedness principle 
jjro(t)jj = as t ^ 00 . This implies that a;o(To) < 0, which is 

a contradiction. Hence there is no hope of finding a rate of energy decay 
which is valid for all x G X] see [141 Lemma 3.1.7] for a related result. 
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(c) The proof of [ini Theorem 2.4] can easily be extended, simply by nsing 
the semigronp property, to show that condition (i) in Theorem 14.11 is 
equivalent to having ||r(t)a:|| = as t —)• oo for all x G D{A^) 

and all integers k > 1. Thus the proof of Theorem 14.21 shows that 
Ex{t) = as t —)• oo for all x G D{A^), k > 1, which is to say that 

smoother orbits have faster energy decay. 


5. The case of Dirichlet boundary conditions 


As mentioned in Section [H the above approach can be modified straight¬ 
forwardly to deal with the model 


wtidt) = 

= 0 , 

ut{0,t) = w{0,t), 
u{C,0) = u{C), 
w{^,0) = wd), 


w(l,t) = 0, 
u^(0, t) = wdo, t) 


^G (-1,0), t >0, 
Cg (0,1), t >0, 
t > 0, 
t > 0, 

ee (-1,0), 

ee (0,1), 


where Neumann boundary condition ttg(—l,t) = 0 appearing in (11.11) has 
been replaced by the Dirichlet boundary condition = 0 for all t > 0. 

Indeed, in this case it is easy to show by means of the Lumer-Phillips theorem 
that the operator A with the new domain 

D{A) = {(tt,n,u;) G Y : u(—1) = n(—1) = u;(l) = 0, 

n(0) = w{0), u'{0) = rc'(O)} 


generates a contraction semigroup on the Hilbert space 
Z = {(u,u,tc) G X : u(—1) = O}. 

Moreover, the spectrum of A again consists only of isolated eigenvalues of 
finite multiplicity, and in fact 

a{A) = |A G C\{0} : \/Asinh(A) cosh(\/A) + cosh(A) sinh(\/A) = O}. 

In particular (t(A) C C_. Arguments completely analogous to those pre¬ 
sented in Section [3] again lead to the optimal resolvent bound ||i?(fs, A)|| = 
0(|s|l/2) 

as jsj —)• oo. This time Theorem O can be applied directly to 
{T{t))t>o, and the result Ex{t) = o(t“^) as t —)• oo for all x G D{A) follows 
at once from the fact that, by the Poincare inequality, the energy defines an 
equivalent norm on Z. 

These results are obtained in [TB] for a similar problem involving an al¬ 
ternative coupling condition and the corresponding more restrictive choice 
of D{A). The technique used in [16] relies on the theory of Riesz spectral 
operators and requires a very detailed spectral analysis. The present ap¬ 
proach based on Theorem 14.11 is rather simpler and in particular requires 
no particular knowledge about those eigenvalues which do not approach the 
imaginary axis, nor any knowledge about any eigenfunctions. 



14 


CHARLES BATTY, LASSI PAUNONEN, AND DAVID SEIFERT 


References 

[1] W. Arendt, C.J.K. Batty, M. Hieber, and F. Neubrander. Vector-valued Laplace trans¬ 
forms and Cauchy problems. Birkhauser, Basel, second edition, 2011. 

[2] G. Avalos and R. Triggiani. The coupled PDE system arising in fluid/structure in¬ 
teraction. I. Explicit semigroup generator and its spectral properties. In Fluids and 
waves, volume 440 of Contemp. Math., pages 15-54. Amer. Math. Soc., Providence, 
RI, 2007. 

[3] G. Avalos and R. Triggiani. Mathematical analysis of PDE systems which govern 
fluid-structure interactive phenomena. Bol. Soc. Paran. Mat., 25:17-36, 2007. 

[4] G. Avalos and R. Triggiani. Semigroup well-posedness in the energy space of a 
parabolic-hyperbolic coupled Stokes-Lame PDE system of fluid-structure interaction. 
Discr. Cont. Dynam. Sys., 2:417-447, 2009. 

[5] G. Avalos. The strong stability and instability of a fluid-structure semigroup. Appl. 
Math. Optim., 55 (2): 163-184, 2007. 

[6] G. Avalos and M. Dvorak. A new maximality argument for a coupled fluid-structure 
interaction, with implications for a divergence-free finite element method. Appl. Math. 
(Warsaw), 35(3):259-280, 2008. 

[7] G. Avalos and R. Triggiani. Rational decay rates for a PDE heat-structure interaction: 
a frequency domain approach. Evol. Equ. Control Theory, 2(2):233-253, 2013. 

[8] G.J.K. Batty, R. Ghill, and Y. Tomilov. Fine scales of decay of operator semigroups. 
J. Europ. Math. Soc., to appear (http://arxiv.org/abs/1305.5365). 

[9] G.J.K. Batty and T. Duyckaerts. Non-uniform stability for bounded semi-groups in 
Banach spaces. J. Evol. Equ., 8:765-780, 2008. 

[10] A. Borichev and Y. Tomilov. Optimal polynomial decay of functions and operator 
semigroups. Math. Ann., 347:455-478, 2010. 

[11] W. Desch, E. Fausangova, J. Milota, and G. Propst. Spectrum of a viscoelastic bound¬ 
ary damping model. J. Integral Equations Appl., 23(4):521-539, 2011. 

[12] T. Duyckaerts. Optimal decay rates of the energy of a hyperbolic-parabolic system 
coupled by an interface. Asymptot. Anal., 51(l):17-45, 2007. 

[13] J. Rauch, Xu Zhang, and E. Zuazua. Polynomial decay for a hyperbolic-parabolic 
coupled system. J. Math. Pures Appl, 84(4):407 - 470, 2005. 

[14] J.M.A.M. van Neerven. The Asymptotic Behaviour of Semigroups of Linear Opera¬ 
tors. Birkhauser, Basel, 1996. 

[15] J.L. Vazquez and E. Zuazua. Large time behavior for a simplified ID model of fluid- 
solid interaction. Comm. Partial Differential Equations, 28(9-10):1705-1738, 2003. 

[16] Xu Zhang and E. Zuazua. Polynomial decay and control of a 1-d hyperbolic-parabolic 
coupled system. J. Differential Equations, 204(2):380-438, 2004. 

[17] Xu Zhang and E. Zuazua. Long-time behavior of a coupled heat-wave system arising 
in fluid-structure interaction. Arch. Ration. Mech. Anal., 184(1):49-120, 2007. 

St John’s College, St Giles, Oxford 0X1 3JP, United Kingdom 
E-mail address: charles.batty@sjc.ox.ac.uk 

Department of Mathematics, Tampere University of Technology, P.O. Box 

553, 33101 Tampere, Finland 

E-mail address: lassi.paunonen@tut.fi 

St John’s College, St Giles, Oxford 0X1 3JP, United Kingdom 
E-mail address: david.seifert@sjc.ox.ac.uk 



